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Abstract 

We  consider  a  thick  porous  metal-hydride  electrode  in  a  temporarily  interrupted  discharge.  During  a  break  the  hy  drogen  and  electrolyte 
concentrations  tend  to  level  steady  states,  which  leads  to  a  recovery  phenomenon  seen  in  the  decrease  in  potential  loss  in  the  electrode  when 
the  discharge  is  restarted.  A  mathematical  model  is  used  in  studying  the  electrochemical  system.  The  potential  loss  is  expressed  as  a  sum  of 
potential  losses  due  to  various  sources,  and  spatial  averaging  is  used  in  analyzing  the  individual  losses  as  well  as  die  reaction  conditions. 
Furthermore,  a  simple  model  is  developed  for  describing  the  behavior  of  the  hydrogen  concentration  during  a  break. 
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1.  Introduction 

We  are  working  on  the  development  of  a  battery  with  a 
metal-hydride  (MH)  electrode  for  electric-vehicle  use.  In  this 
kind  of  application  the  battery  is  discharged  with  a  strongly 
varying  current.  Furthermore,  the  discharge  is  frequently 
completely  interrupted.  In  this  paper  we  investigate  the  per¬ 
formance  of  a  thick  MH  electrode  over  a  discharge  break 
using  a  mathematical  model. 

Mathematical  models  for  MH  electrodes  are  presented  in 
Refs.  [  1-3],  for  example.  A  complete  cell  with  a  thin  MH 
electrode  is  modeled  in  Ref.  [4].  Recovery  or  relaxation 
phenomena  for  cells  with  different  electrodes  have  been 
considered  in  Refs.  [5,6] . 

In  this  paper,  we  use  the  model  developed  in  Ref.  [  3  ] ,  but 
reconsider  the  potential  losses  in  more  detail,  and  keeping  in 
mind  the  interrupted  discharge.  As  before,  we  use  spatial 
averaging  to  analyze  the  potential  losses  due  to  various  mech¬ 
anisms.  Moreover,  the  averaging  method  is  applied  in  stud¬ 
ying  the  reaction  conditions,  such  as  the  electrolyte  and 
hydrogen  concentrations  and  the  reaction  rate  density. 

Our  mathematical  model  is  simple  so  that  analytic  methods 
can  be  applied  more  easily  in  investigating  the  simulations 
results.  On  the  other  hand,  with  a  simple  model  one  looses 
some  accuracy;  this  model  should  only  be  used  to  study  and 
interpret  phenomena  and  their  relationship  in  a  thick  MH 
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electrode  instead  of  trying  to  produce  pr  ecise  predictions  for 
specific  quantities. 

2.  Mathematical  model 

A  brief  description  of  the  mathemaucal  model  is  given; 
more  details  can  be  found  in  Ref.  [3] . 

2.1 .  Porous  MH  electrode  and  the  cell 

The  porous  MH  electrode  is  manufactured  by  bonding  MH 
particles  together  with  polytetrafluoroeihylene  (PTFE)  and 
carbon.  In  the  cell  the  MH  electrode  is  immersed  in  KOH 
electrolyte. 

The  void  volume  fraction,  i.e.  the  porosity  of  the  electrode 
is  denoted  by  e.  Let  v  be  the  volume  fraction  of  the  MH  alloy 
in  the  electrode.  Due  to  the  passive  bonding  material 
e  +  v  <  1 .  The  porosity  and  tortuosity  of  the  electrode  are 
taken  into  account  by  multiplying  the  ionic  diffusion  coeffi¬ 
cients  Dfa*  and  conductivities  cr^  by  el  5 .  The  subscript 
‘free’  refers  to  solution  outside  any  porous  structure.  We 
define  the  effective  ionic  diffusion  coefficients  and  conduc¬ 
tivities  by 

D±=e<  5Df*e,  <J±  =  e‘  Vf±e  ( 1 ) 

In  this  article,  we  concentrate  on  the  MH  electrode. 
Accordingly,  the  counter  electrode  is  viewed  only  as  a  source 
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Current  Collector  Mesh 


0  x,  xc  xr 

Fig.  1 .  Schematic  diagram  of  the  electrode  system.  The  section  of  the  cell 
tube  with  the  electrodes  is  shown  straightened  here.  The  circular  tube  is 
formed  when  the  point  x  =  L  (far  right,  not  shown)  is  connected  with  the 
point  x  —  0. 

of  ions  and  modeled  as  very  thin.  In  laboratory  experiments 
a  nickel  mesh,  for  example,  can  be  the  thin  counter  electrode. 
Its  location  is  denoted  by  x{.  The  locations  of  the  left  and  right 
ends  of  the  MH  electrode  are  denoted  by  x}  and  xr,  respec¬ 
tively.  The  thin  current  collector  is  placed  at  the  mid  point  of 
the  MH  electrode  at  xc  =  (jq + xT)  / 2,  see  Fig.  1 . 

We  use  the  same  cell  geometry  as  in  Refs.  [3,7] ,  that  is, 
we  consider  a  circular  cell.  Note  that  this  is  only  a  mathe¬ 
matical  model  and  not  a  description  of  a  real  practical  cell. 
As  we  make  L,  the  total  length  of  the  cell  tube,  large  compared 
with  xT  —  xh  the  thickness  of  the  MH  electrode,  the  situation 
corresponds  to  two  electrodes  in  a  large  reservoir  of  electro¬ 
lyte.  This  is  a  typical  experimental  setup  in  testing  MH  elec¬ 
trodes.  To  avoid  discontinuities  in  the  diffusion  coefficients 
and  the  conductivities  we  assume  that  outside  the  electrodes 
the  tube  is  filled  with  a  passive  porous  material  having  the 
same  porosity  as  the  MH  electrode. 

Finally,  the  discharge  flux  density  is  denoted  by  J. 

2.2.  Reaction  kinetics 


The  operation  of  an  MH  electrode  is  based  on  the  reaction 
MHads  +  OH ~  +H20  +  M  (2) 


We  use  a  very  simple  Butler-Volmer-type  expression  [8] 
to  describe  the  reaction  kinetics  of  Eq.  (2).  Let  c  and  hs  be 
the  electrolyte  concentration  and  the  concentration  of  hydro¬ 
gen  on  the  surface  of  a  MH  particle,  respectively.  The  con¬ 
centration  of  water  is  assumed  to  be  constant  and  denoted  by 
cH20.  The  reaction  rate  density,  r  (mol/  (m3  s) ),  can  then  be 
modeled  by 


r(t,x)  = 


P  exp*^ - — - •|c(U)/is(r,x) 

f  0.5  FV(tyx)\  1 

-  kc  exp« - — - >cH2 o  J 


V(tj)  =  (p(t,X )  -  (^(tyX )  (4) 

is  the  difference  between  the  potential  in  the  solid  and  the 
potential  in  the  electrolyte,  denoted  by  cp  and  </>,  respectively, 
measured  with  a  Hg/HgO  reference  electrode.  Furthermore, 
ka  and  kc  are  the  reaction  rate  constants  and  T is  temperature, 
which  is  assumed  to  be  constant.  Since  we  measure  exchange 
current  density  in  A/g  we  have  included  the  factor  vplF , 
where  p  is  the  mass  density  of  the  MH  alloy,  in  Eq.  (3). 


2.5.  Ion  transport  in  the  electrolyte 


The  ions  are  driven  by  diffusion  and  migration  in  the  elec¬ 
trolyte.  The  motion  of  the  ions  is  modeled  with  a  dilute  solu¬ 
tion  theory  from  Ref.  [8]  using  constant  ionic  diffusion 
coefficients  and  ionic  conductivities.  The  electrolyte  concen¬ 
tration  c  and  potential  4>  satisfy  the  equations 


a  ( 

€Cf(tyX)  =  ““I 


D  cx(tyx)  + 


ect(tyx)  = 


-s(-° 


) 


+  £(*,*) 


cx(t,x)  — =r  </>,(*>*) 


(5) 

(6) 


which  describe  the  motion  of  the  OH'  and  K+  ions, 
respectively. 

The  inhomogeneous  term,  £,  in  Eq.  (5)  represents  the 
effect  of  the  electrode  reactions.  OH  ions  are  produced  on 
the  thin  counter  electrode  with  the  rate  corresponding  to  the 
discharge  flux  density,  J ,  and  consumed  in  the  MH  electrode 
with  the  rate,  r,  given  by  Eq.  (3).  Accordingly 


Eit^^Jm^W-ritpc)  (7) 

where  a  Dirac  delta  function  8^  is  used  in  describing  the  thin 
counter  electrode.  The  K+  ions  do  not  take  part  in  the  reac¬ 
tions  and,  hence,  there  is  no  corresponding  term  in  Eq.  (6) . 
Due  to  the  conservation  of  charge  the  condition 


jr(t&de=m  (8) 


must  hold  for  all  t  >  0. 

Mathematically  the  circular  tube  is  described  by  periodic 
boundary  conditions.  We  require  that  the  electrolyte  concen¬ 
tration  and  potential  and  the  electric  field  at  x  —  L  coincide 
with  the  their  respective  values  at  x  -  0.  More  precisely 

c(f,0)  =c(t,L)  (9) 

</>(r,0)  =  4>(f,L)  (10) 

(k(r,0  )=cf>x(t,L)  (ID 

We  emphasize  that  these  quantities  are  not  periodic  over  the 
MH  electrode  but  over  the  whole  cell  tube. 

In  the  beginning  of  the  discharge  the  electrolyte  concen¬ 
tration  is  constant  everywhere.  Therefore,  we  have  the  initial 
condition 


where 


c(0,x)  =c0. 


(12) 
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2.4.  Hydrogen  diffusion  in  MH  material 


and 


The  hydrogen  diffusion  is  modeled  in  a  simplified  geom¬ 
etry.  By  scaling  the  diffusion  coefficient  we  make  sure  that 
the  time  scale  of  the  hydrogen  diffusion  remains  correct. 

Let  Z,  the  maximal  diffusion  distance  in  the  model,  be 

Z=  v!  y  (13) 

where  y  is  the  active  surface  area  density  (m2/m3)  of  the 
MH  alloy.  Next,  we  define  the  scaled  diffusion  coefficient,  d , 
by 

d=Tff!  td jff  —  v 2 1  ( y^diff)  (14) 


<px{tJr)  =0  (23) 

state  that  the  ends  of  the  electrode  are  insulated. 

3.  Potential  losses  and  reaction  conditions 

We  define  the  electrode  potential  <P  to  be  the  potential 
difference  between  the  current  collector  and  the  electrolyte 
between  the  counter  and  MH  electrodes 

<P(t)  =  <p(tyX c)  -  <£(L*i)  (24) 


where  t^ ff  is  the  characteristic  time  of  the  hydrogen  diffusion. 
Finally,  the  characteristic  time  of  the  hydrogen  diffusion  is 
defined  by 

(15) 

where  rMH  is  a  typical  diffusion  distance  in  real  MH  particles 
and  Dh  is  the  diffusion  coefficient. 

Designate  the  concentration  of  hydrogen  in  the  MH  mate¬ 
rial  by  h .  We  model  the  hydrogen  diffusion  with  the  basic 
diffusion  equation 

ht(tyzj)  =dhzz(ttz*x)  ( 16) 

where  t>0,0<z<Z  and  In  this  section  x  is  only 

a  parameter.  The  boundary  condition 

ydhz(t,Z,x)  =  -  r(t,x)  ( 17) 

models  the  hydrogen  flow  through  surfaces  of  particles  due 
to  the  reaction  (2).  The  other  boundary  condition 

dhz(ty  0,x)=0  (18) 


The  electrode  potential  is  measured  with  a  Hg/HgO  ref¬ 
erence  electrode.  The  local  momentary  equilibrium  potential 
V0  and  the  steady-state  open-circuit  equi  1  ibrium  potential  V0>0O 
are  obtained  from  Eq.  (3)  by  setting  r  ^0  and  using  appro¬ 
priate  electrolyte  and  hydrogen  concentrations 


ch2o 

c(tj)hs(tj) 


(25) 


and 


H20 


Co(l-f/fmax)k 


■} 


(26) 


where 


'max 


v(xr-xx)h0 

J 


(27) 


is  the  maximal  theoretical  discharge  time  when  the  electrode 
is  discharged  with  the  constant  flux  density,  J.  If  the  current 
is  switched  off  at  t  =  tc  then 


corresponds  to  a  symmetry  condition  in  the  middle  of  a  real 
particle.  We  assume  that  the  initial  hydrogen  concentration 
is  constant  throughout  the  MH  electrode 


limVo(Lx)  =  V0>00(rc)  (28) 

t  — ►  00 

We  may  define  the  total  potential  loss  in  the  electrode 


h(0,zj)=ho  (19) 

Finally,  the  surface  concentration  of  hydrogen,  /is,  is 
defined  by 

K{tjc)  =h(t,Z,x)  (20) 

2.5.  Electric  current  in  the  electrode  material 

Using  Ohm’s  law  we  obtain  an  equation  for  the  potential 
<p  in  the  electrode  material 

K(pxx(t,x)  =  F(r(tyx)  -J(t)hXc(x))  (21) 

where  k  is  the  effective  electric  conductivity  of  the  electrode 
material.  The  thin  current  collector  at  xc  =  (xt + xr)  /2  is  mod¬ 
eled  with  a  Dirac  delta  function  hXc.  In  Eq.  (21)  JCj<x<xr 
and  t  is  only  a  parameter. 

The  boundary  conditions 

(pAtyXi)  =0  (22) 


vloss(0  =  ^(0-v0iOC(0  (29) 

As  in  Ref.  [3]  we  divide  the  total  potential  loss  into  four 
components:  (i)  the  potential  loss  in  the  electrolyte,  rjj;  (ii) 
the  reaction  overpotential,  tj;  (iii)  the  potential  loss  in  the 
solid  phase  tjs,  and  (iv)  the  loss  due  to  concentration  polar¬ 
ization  (the  deviation  of  the  local  momentary  equilibrium 
potential  from  the  steady-state  open-circuit  equilibrium 
potential),  7jp.  Both  electrolyte  and  hydrogen  concentration 
polarization  contribute  to  7}p.  To  be  more  precise 


i?i(M 0  =  4>{t>x)  -  <f>(tjcx)  (30) 

7 i(tjc)=V(tj)-V0(t,x)  (31) 

y h(tpc)  =  <p(tjcG )  -  <p(tjc)  (32) 

and 

V p(t,x)  =  Vo (tyX)  -  V0,oo(r)  (33) 


To  describe  the  total  significance  of  a  individual  loss  com¬ 
ponent  at  a  given  time  with  a  single  number  we  first  compute 
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the  total  power  loss  per  projected  area  (W/m2)  due  to  the 
loss  source  in  question.  For  example,  consider  the  potential 
loss  in  the  electrolyte.  Since  the  total  length  of  the  cell  tube 
is  large  compared  with  the  thickness  of  the  MH  electrode,  it 
is  reasonable  to  assume  here  that  the  current  outside  the  MH 
electrode  is  negligible.  The  current  density  for  xx  <  x  <  xT  is 

X 

Jy-  Jk  t,f))  (34) 

X I 


so  that  the  total  power  loss  in  the  MH  electrode  due  to  the 
resistance  of  the  electrolyte  is 

Xr  /  *  \ 

xi  xt 

(  t  \  Xr 

=  v  ~  J r(f,Tj)dTjp(^>  jFr^€) 

XI  XI 


XT 

=  (Fr(t,&<Kt,£)dZ-FJ<Htxi) 

X] 


Xr 

XI 


Xr 

=  jFrit&rhit&di-  (35) 

XI 

Now  the  corresponding  average  potential  loss  in  the  electro¬ 
lyte  %  is  obtained  by  dividing  the  power  loss  density  by  the 
current  density 

Xr  Xr 

jFr(tSi 

TO=a~s - - - - -  (36) 

JVr(/,£)  df 

X\ 

The  remaining  three  average  potential  losses  %  7^,  and 
are  computed  similarly.  We  use  the  above-defined  average 
quantities  to  explain  the  changes  in  the  total  potential  loss  in 
the  MH  electrode  over  a  break  in  discharge. 

We  also  compute  the  corresponding  weighted  averages  of 
the  electrolyte  concentration  and  the  surface  concentration  of 
hydrogen,  denoted  by  c  and  JTS: 


Jr(r,£)c(f,£)  df 

cU)=- - - -  (37) 


Xr 

jV(r,£)/j,(f,£)d£ 

K(D=- - - -  (38) 

Although  there  are  no  direct  physical  interpretations  in  these 
cases,  these  averaged  quantities  give  us  a  clue  to  the  time 
evolution  of  ‘effective’  reaction  conditions.  Furthermore,  we 
define  the  average  reaction  rate  density,  7,  by 

Xt 

jlr(t,&}2  df 

Ht)  =~ - - -  (39) 


This  quantity  reflects  the  localization  of  the  reaction.  The 
larger  the  value  of  r  the  more  localized  the  reaction  is.  It  is 
easy  to  show  that  r(t)  >J/  (*r  —  x})  for  all  f,  where  the  lower 
limit  corresponds  to  the  situation  in  which  the  reaction  rate 
density  is  constant  throughout  the  MH  electrode,  and  the 
reaction  is  fully  spread  out. 

In  addition  to  the  fact  that  these  averaged  quantities 
describing  the  reaction  conditions  are  interesting  as  such  they 
can  be  used  to  study  the  averaged  potential  losses  further.  We 
apply  the  average  electrolyte  and  hydrogen  concentrations 
and  the  average  reaction  rate  density  in  analyzing  the  average 
reaction  overpotential  and  concentration  of  polarization  loss. 
First  we  define  the  reference  values  for  the  reaction  rate  den¬ 
sity,  the  surface  concentration  of  hydrogen  and  the  electrolyte 
concentration. 

Let  the  reference  reaction  rate  density,  r0,  be 

>o  =  — ^ —  (40) 

xr-*i 

That  is,  r0  corresponds  to  the  situation  where  the  reaction  rate 
density  is  constant  throughout  the  MH  electrode.  Let  h^{t) 
be  the  concentration  that  results  if  the  total  amount  of  hydro¬ 
gen  in  the  electrode  at  time  t  is  evenly  spread  out.  In  other 
words 


hao(t)  =  (l-t/tmax)h0 


(41) 


The  reference  value  for  the  electrolyte  t  oncentration  is  chosen 
to  be  c0.  In  the  following  the  electrolyte  concentration,  the 
surface  concentration  of  hydrogen  and  the  reaction  rate  den¬ 
sity  are  considered  relative  to  their  reference  values  c0,  hoc 
and  r0,  respectively. 

We  start  with  the  reaction  overpotcntial.  Assume  that  the 
cathodic  current  is  negligible  during  the  discharge,  we  then 
easily  obtain  from  Eq.  (3) 


V(^)- 


RT 


'44 

l¥a 


r(t,x) 


0.5  F  {vpkac(tjc)hs(t^ c) 


} 


(42) 


Combining  Eq.  (42)  with  Eq.  (25)  yields  after  some  manip¬ 
ulations  an  approximation  for  the  reaction  overpotential 


V(tjc)  ~ 


RT 

F 


F2  (r(M  ')2 
v2p2kJiccH2Qc{tj,)h,  tjc) 


and 


(43) 
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By  linearizing  the  previous  expression  with  respect  to  the 
triple  (r0icQ,hoo(t))  and  substituting  the  result  into  the 
definition  of  the  average  loss  we  get 


where  we  have  ignored  the  residual  term. 

Similarly,  the  concentration  polarization  loss  can  be 
linearized 


(45) 


Now  we  can  use  the  averaged  concentrations  and  the  reac¬ 
tion  rate  density  in  Eqs.  (44)  and  (45)  in  explaining  the 
behavior  of  the  reaction  overpotential  and  the  concentration 
polarization  loss. 

By  using  the  averaging  technique  and  removing  the  spatial 
dimension  we,  in  a  way,  approximate  the  thick  MH  electrode 
with  an  infinitely  thin  one  in  which  the  averaged  concentra¬ 
tions  and  the  reaction  rate  density  are  the  prevailing  reaction 
conditions. 

We  emphasize  that  the  approximative  formulas  derived 
above  are  presented  only  to  help  us  to  interpret  the  simulation 
results. 


Table  1 

General  parameters 


Parameter 

Value 

References 

DP 

10"'2  m2/s 

m 

f^firee 

3.33  X  10"9m2/s 

[9] 

At use 

1.07  X  10-2  S  m2/mol 

[9] 

t  + 

0.22 

[10] 

2.13X  10-*  m2/s 

derived 

f^free 

7.57  X10~5  m2/s 

derived 

14.12  S/m 

derived 

50.08  S/m 

derived 

K 

4.4  X  10 ~4  m6/(kg  mol  s) 

derived 

kc 

5.7 X  10" 16  m3/ (kg  s) 

derived 

K 

150  S/m 

measured 

C0 

6000  mol/m3 

CH20 

40000  mol/m3 

ho 

77400  mol/m3 

J 

0.01  mol/(m2  s) 

L 

0.1  m 

xt 

0.001  m 

Xr 

0.005  m 

y 

8  X  104  m2/m3 

V 

0.4 

e 

0.5 

P 

7.74  X103  kg/m3 

r0 

2.5  mol/ (m3  s) 

J/(xr-^) 

rMH 

15  X  10-6  m 

T 

293  K 

rate.  The  reference  reaction  rate  density  is  then  2.5  mol/ 
(m3  s). 


4.  Simulation 


4.2.  Results 


Here  we  consider  a  case  where  the  discharge  is  temporarily 
interrupted.  The  simulations  predict  a  recovery  phenomenon 
in  the  potential  loss. 

4.1.  Parameters 

The  values  of  the  parameters  are  given  in  Table  1.  The 
concentration  of  water  and  the  ionic  diffusion  coefficients 
and  conductivities  are  determined  for  6  mol/1  KOH  electro¬ 
lyte.  The  reaction  rate  constants  are  computed  assuming  a 
fully  charged  electrode  and  the  initial  electrolyte  concentra¬ 
tion  together  with  a  measured  equilibrium  potential  value  of 
—  0.927  V  versus  Hg/HgO  and  a  measured  exchange  current 
density  of  210  mA/g  from  Ref.  [11]. 

The  initial  concentration  of  hydrogen  is  1  wt.%  in  the  MH 
alloy  LMNi3 6Co0 7A10 4Mno 3  (Treibacher,  Austria).  The 
porosity  and  the  MH  material  volume  fraction  as  well  as  the 
thickness  of  the  electrode  are  known  from  the  preparation  of 
prototype  electrodes.  The  surface  area  density  and  the  typical 
diffusion  distance  correspond  to  a  precycled  electrode. 

The  total  length  of  the  cell  tube,  L,  is  chosen  large  com¬ 
pared  with  xr — xh  the  thickness  of  the  MH  electrode,  so  that 
the  current  outside  the  electrodes  is  negligible. 

The  electrode  is  discharged  with  a  discharge  flux  density 
of  0.01  mol/m2  s  (100  mA/cm2)  corresponding  to  C/3.3 


In  this  simulation  the  discharge  is  interrupted  at  6000  s 
(the  corresponding  state-of-discharge  is  about  50%)  and  con¬ 
tinued  after  a  break  of  7200  s. 

Fig.  2  depicts  Vloss,  the  total  potential  loss  in  the  electrode. 
It  can  be  seen  that  the  potential  loss  immediately  after  the 
break  is  smaller  than  immediately  before  it.  The  difference  is 
about  10  mV. 


Fig.  2.  Total  potential  loss  in  the  MH  electrode,  Vloss.  The  curves  before  and 
after  the  break  are  separated  with  a  dashed  line 
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Fig.  3.  Averaged  potential  losses  ^,1 ?s,i7p,  and  17.  The  curves  before  and  after 
the  break  are  separated  with  a  dashed  line. 

In  this  case,  the  recovery  is  easily  explained  with  Fig.  3 
showing  the  averaged  potential  losses.  The  most  important 
contribution  to  the  decreased  loss  is  the  potential  loss  in  the 
electrolyte,  which  decreases  about  10  mV.  The  changes  in 
reaction  overpotential  and  concentration  polarization  loss 
cancel  each  other  and  the  increase  in  potential  loss  in  the  solid 
is  negligible. 

We  take  now  a  closer  look  in  the  individual  loss  compo¬ 
nents.  During  the  break  the  hydrogen  concentration  con¬ 
verges  towards  a  level  value,  see  Fig.  4.  At  the  same  time  the 
electrolyte  concentration  also  levels  out  as  shown  in  Fig.  5. 
Thus  the  situation  after  the  break  resembles  the  start-of-dis- 
charge  and  the  reaction  zone  starts  again  from  the  counter 
electrode,  Fig.  6. 

The  decrease  in  potential  loss  in  the  electrolyte  is  caused 
by  the  reaction  zone  moving  back  to  the  counter  electrode 
and  the  resulting  shortening  of  the  path  of  ions  in  the  electro¬ 
lyte.  The  potential  loss  in  the  solid  increases  since  the  distance 
between  the  reaction  zone  and  the  current  collector  gets 
longer. 

We  use  the  linear  approximations  (44)  and  (45)  and  the 
relative  average  electrolyte  concentration  c/c0,  hydrogen 
concentration  hjh and  the  relative  average  reaction  rate 


Fig.  4.  Relative  surface  concentration  of  hydrogen,  hjh0.  During  the  break 
the  concentration  tends  to  a  level  value  of  about  0.5. 


Fig.  5.  Relative  electrolyte  concentration,  c/< ,,  During  the  break  the  con¬ 
centration  tends  to  the  level  initial  value  c0. 


Fig.  6.  Relative  reaction  rate  density,  r/ r0.  Thv  negative  values  during  the 
break  indicate  charging  related  to  the  leveling  of  hydrogen. 

density  ft  r0  shown  in  Fig.  7  to  analyze  the  changes  in  reaction 
overpotential  and  concentration  polarization  loss.  Generally, 
the  reaction  rate  density  varies  the  most  while  the  changes  in 
electrolyte  and  hydrogen  concentrations  are  considerably 
smaller. 

In  the  beginning  of  the  discharge  the  reaction  almost  com¬ 
pletely  takes  place  next  to  the  counter  electrode.  This  can  be 
seen  in  the  large  average  reaction  rate  density.  With  the  time 
the  reaction  spreads  deeper  into  the  electrode  and  at  t— 6000  s 


Fig.  7.  Relative  values  of  averaged  reaction  rate  density,  r/r0,  electrolyte 
concentration,  c/c0,  and  surface  concentration  of  hydrogen,  hJK 0.  The 
curves  before  and  after  the  break  are  separated  with  a  dashed  line. 
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the  reaction  rate  density  is  almost  uniform  as  indicated  by  the 
value  of  about  1  of  the  relative  average  reaction  rate  density. 
During  the  discharge,  the  electrolyte  concentration  decreases 
away  from  the  counter  electrode.  Thus,  as  the  reaction  moves 
further  into  the  electrode  it  experiences  lower  electrolyte  con¬ 
centration  as  seen  in  Fig.  7.  In  the  relative  average  hydrogen 
concentration,  we  see  the  effect  of  the  decreasing  total  supply 
of  hydrogen  causing  the  increase  observed  before  the  break. 

From  Figs.  3,  6,  and  7  and  Eq.  (44)  it  is  easy  to  see  that 
the  reaction  rate  density  dominates  the  behavior  of  the  reac¬ 
tion  overpotential.  Before  the  break,  the  reaction  takes  place 
in  the  middle  of  the  MH  electrode  and  is  almost  completely 
spread  out.  At  the  restart  of  the  discharge  the  reaction  is  again 
relatively  localized  next  to  the  counter  electrode  and  an 
increase  in  reaction  overpotential  is  developed.  The  increase 
in  electrolyte  concentration  and  the  decrease  in  surface  con¬ 
centration  of  hydrogen  have  a  decreasing  and  an  increasing 
effect  on  the  reaction  overpotential,  respectively. 

In  the  concentration  polarization  loss  there  are  two  com¬ 
peting  factors  and  the  change  is  small.  It  appears  that  here  the 
increase  in  electrolyte  concentration  overrides  the  effect  of 
the  decreasing  surface  concentration  of  hydrogen.  As  a  result, 
a  3  mV  decrease  in  the  concentration  polarization  loss  is 
produced. 

4.3.  Hydrogen  concentration  during  a  break 


Fig.  8.  Comparison  between  the  simulated  and  the  approximated  surface 
concentration  of  hydrogen  (Eq.  (48 ) )  next  to  the  c ounter  electrode  at  x = 
during  the  discharge  break.  The  dashed  vertical  line  marks  the  time  when 
the  discharge  was  interrupted  and  the  dashed  horizontal  line  corresponds  the 
hydrogen  concentration  with  the  remaining  hydrogen  at  f=6000  s  evenly 
distributed  in  the  electrode. 


(47) 


It  is  easy  to  see  that  the  solution  is 


K (t>x)  =M*c)  +  (hs(tcyx)  -  MO)  ex 


(48) 


We  now  study  closer  the  most  important  factor  in  the 
electrode  recovery,  the  hydrogen  concentration  change  dur¬ 
ing  a  discharge  break.  The  leveling  of  the  hydrogen  shown 
in  Fig.  4  allows  the  reaction  to  start  again  from  the  counter 
electrode,  which  leads  to  the  decrease  in  total  potential  loss. 

Let  tc  be  the  time  when  the  discharge  is  interrupted.  The 
corresponding  level  steady- state  hydrogen  concentration 
hoo(tc)  is  defined  by  Eq.  (41).  That  is,  /ioo(fc)  is  the  concen¬ 
tration  towards  which  the  hydrogen  concentration  in  the  elec¬ 
trode  converges  after  the  discharge  has  been  interrupted  at 
time  tc. 

To  estimate  roughly  the  rate  of  convergence  of  the  hydro¬ 
gen  we  assume  that  the  diffusion  and  the  migration  in  the 
electrolyte  are  infinitely  fast  and  there  is  no  electric  resistance 
in  the  solid.  We  also  assume  that  the  hydrogen  diffusion  in 
the  MH  particles  is  infinitively  fast  and  that  the  electrolyte 
concentration  and  the  potential  difference  between  the  solid 
and  the  electrolyte  take  their  steady-state  values  c0  and  V0.°° 
immediately  after  the  current  has  been  switched  off.  We  also 
assume  that  the  hydrogen  diffusion  in  the  MH  particles  is 
infinitely  fast.  With  these  assumptions  the  surface  concentra¬ 
tion  of  hydrogen  satisfies  for  t  >  tc 

a,  .  p|~,  fo.5FV0M\  .  _ 

=  --|  K exP"| - — - j-CfAM 

f-0.5FVoM\  1 

-  kc  exp< - — - >cHl0  I  ( 46 ) 


where  the  characteristic  time,  r,  is  given  by 


(49) 


Hence  the  hydrogen  concentration  converges  exponentially 
towards  the  steady  state  h 

Substituting  the  values  for  the  parameters  we  find  that  if 
fc  =  6000  s  corresponding  to  the  state  of  discharge  of  about 
50%  then  r*  3500  s,  which  fits  the  results  reasonably  well. 
A  comparison  of  the  approximation  and  the  simulated  surface 
concentration  of  hydrogen  is  presented  in  Fig.  8,  in  which  the 
situation  right  next  to  the  counter  electrode  is  investigated. 
We  notice  that  the  characteristic  time  given  above  appears  to 
be  too  short,  that  is,  the  concentration  given  by  Eq.  (48) 
converges  too  rapidly. 

The  main  reason  for  this  is  that  after  the  current  has  been 
cut  there  still  exist  non  zero  electrolyte  concentration  and 
potential  gradients  that  sustain  the  ion  flux  of  the  leveling 
process.  As  a  result  the  electrolyte  concentration  and  potential 
differ  from  their  level  steady-state  values.  Taking  this  into 
account  in  Eq.  (46)  and  solving  for  the  surface  concentration 
of  hydrogen  corresponding  to  the  balance  condition 


d 

T*s(U)  =0 
or 


(50) 


we  obtain 


CClhO 

Kc{tjc) 


FV(t,x)\ 
RT  ) 


where 


(51) 
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This  hydrogen  concentration  at  a  given  time,  ty  deviates  from 
hco  towards  the  concentration  before  the  break  and  the  devi¬ 
ation  is  naturally  at  its  largest  values  immediately  after  the 
discharge  has  been  interrupted.  Consequently,  the  simulated 
convergence  is  slower  than  the  approximated  one.  Moreover, 
the  lower  the  porosity  of  the  MH  electrode  the  larger  are  the 
electrolyte  concentration  and  potential  gradients.  Thus,  lower 
porosity  leads  to  slower  convergence  of  hydrogen 
concentration. 

Note  that  Tis  a  decreasing  function  of  t  and  that  rit^)  = 
0  s.  This  behavior  is  explained  easily.  The  more  the  electrode 
is  discharged  the  smaller  amount  of  hydrogen  has  to  move  in 
the  process  of  convergence  to  the  constant  steady  state. 

5*  Conclusions 

With  the  average  losses,  see  Ref.  [3] ,  the  average  electro¬ 
lyte  and  hydrogen  concentration,  and  the  average  reaction 
rate  density  presented  in  this  article  it  is  easy  describe  and 
analyze  the  discharge  performance  of  a  thick  MH  electrode. 

These  tools  were  applied  in  studying  a  break  of  discharge. 
It  is  clear  that  the  dominating  factor  in  the  recovery  of  the 
electrode  during  a  break  is  the  decreasing  potential  loss  in  the 
electrolyte. 

The  changes  in  the  degree  of  localization  of  the  reaction 
determine  the  behavior  of  the  second  most  important  source 
of  loss,  the  reaction  overpotential.  In  the  middle  of  the  dis¬ 
charge  the  reaction  is  almost  completely  uniform  but  in  the 
beginning,  and  especially  in  the  end  of  the  discharge,  the 
reaction  is  strongly  localized  causing  an  increase  in  reaction 
overpotential. 

We  developed  a  simple  model  for  estimating  the  leveling 
of  the  hydrogen  concentration  during  a  break,  which  is  the 
basic  phenomenon  in  the  recovery  of  a  thick  MH  electrode. 

6,  List  of  symbols 

c  electrolyte  concentration,  mol  /  m3 
c  average  electrolyte  concentration,  mol  / m3 
c0  initial  concentration  of  the  electrolyte,  mol/ m3 
cH2o  concentration  of  water  in  the  electrolyte,  mol/ m3 
d  pseudo  diffusion  coefficient  in  the  hydrogen 
diffusion  model,  m2/s 

Dfree  integral  diffusion  coefficient  of  KOH,  m2/s 
D+  effective  diffusion  coefficient  of  K  +  ions,  m2/s 
D~  effective  diffusion  coefficient  of  OH-  ions,  m2/s 
Z)free  diffusion  coefficient  of  K+  ions,  m2/s 
Df7ee  diffusion  coefficient  of  OH“  ions,  m2/s 
D H  diffusion  coefficient  of  hydrogen  in  the  MH  alloy, 

m2/s 

F  Faraday’ s  constant,  96  500  C/  mol 

h  hydrogen  concentration  in  MH  alloy,  mol/ m3 
h0  initial  concentration  of  hydrogen  in  MH  alloy, 
mol/m3 


h*.  steady-state  concentration  of  hydrogen,  mol/m3 
hs  surface  concentration  of  hydrogen  in  the  MH  alloy, 
mol/m3 

h&  average  surface  concentration  of  hydrogen  in  the 
MH  alloy,  mol/m3 

J  discharge  flux  density,  mol  / m 2s 

anodic  reaction  rate  constant,  m6/ (kg  mol  s) 
kc  cathodic  reaction  rate  constant ,  m3 /  ( kg  s ) 

L  total  length  of  the  cell,  m 
r  reaction  rate  density  in  the  MH  electrode,  mol/ 

(m3  s) 

7  average  reaction  rate  density  in  the  MH  electrode, 
mol/(m3  s) 

r0  reference  reaction  rate  densit)  in  the  MH  electrode, 
mol/(m3  s) 

rMH  diffusion  distance  in  real  MH  particles,  m 
R  universal  gas  constant,  8.3143  J/ (mol  K) 

t  time,  s 

t +  transference  number  of  K+  ions 

tc  time  when  the  discharge  is  interrupted,  s 
t^ff  characteristic  diffusion  time  in  MH  particles,  s 
tnsn  maximal  discharge  time,  s 

T  temperature,  K 

V  potential  difference  between  the  MH  electrode  and 

the  electrolyte  vs.  Hg/HgO,  V 
Vloss  total  potential  loss  in  the  MH  electrode,  V 
V0  local  momentary  equilibrium  potential  vs.  Hg/HgO, 

V 

Vo oo  steady-state  open-circuit  equil ibrium  potential  vs. 
Hg/HgO,  V 

jc  spatial  coordinate  across  the  e  iectrode,  m 
xc  location  of  the  current  collects  >r,  m 

jc!  location  of  the  counter  electrode  and  the  left  end  of 

the  MH  electrode,  m 

jcr  location  of  the  right  end  of  the  MH  electrode,  m 
z  spatial  coordinate  in  the  hydrogen  diffusion  model, 
m 

Z  maximum  diffusion  distance  in  the  hydrogen 
diffusion  model,  m 

Greek  letters 

y  active  surface  area  density  in  the  MH  electrode, 

1/m 

S€  Dirac’s  delta  function  located  at  x  =  f 
e  porosity  of  the  MH  electrode 

ilfree  equivalent  conductance  of  K(  )H,  S  m2/ mol 
7)  reaction  overpotential,  V 

771  potential  loss  in  the  electrolyte,  V 

7}p  concentration  polarization  loss,  V 

potential  loss  in  the  solid  phase,  V 
I)  average  reaction  o verpotentia  I ,  V 

%  average  potential  loss  in  the  electrolyte,  V 
average  concentration  polarization  loss,  V 
rf$  average  potential  loss  in  the  $<  did  phase,  V 

k  effective  conductivity  of  the  e  Iectrode  material,  S/m 

v  MH  volume  fraction 
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( p  electric  potential  in  the  electrolyte,  V 

<p  electric  potential  in  the  solid  phase,  V 

<P  electrode  potential  vs.  Hg/HgO,  V 
p  density  of  the  MH  material,  kg/ m3 
< effective  ionic  conductivity  of  K+  ions,  S/m 
cr~  effective  ionic  conductivity  of  OH“  ions,  S/m 
crfoe  ionic  conductivity  of  K+  ions,  S/m 
crIVee  ionic  conductivity  of  OH"  ions,  S/m 
r  characteristic  time  of  convergence  of  hydrogen  to 
the  steady  state,  s 
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